Saturation of interband absorption in graphene 
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The transient response of an intrinsic graphene, which is caused by the ultrafast interband tran- 
sitions, is studied theoretically for the range of pumping correspondent to the saturated absorption 
regime. Spectral and temporal dependencies of the photoexcited concentration as well as the trans- 
mission and relitive absotption coefficients are considered for mid-IR and visible (or near-IR) spec- 
tral regions at different durations of pulse and broadening energies. The characteristic intencities 
of saturation are calculated and the results are compared with the experimental data measured for 
the near-IR lasers with a saturable absorber. The negative absorption of a probe radiation during 
cascade emission of optical phonons is obtained. 

PACS numbers: 78.47.jb, 78.67.Wj, 42.65.-k 



I. INTRODUCTION 

The character of nonlinear response under the ultra- 
fast interband excitation of an intrinsic graphene is de- 
termined by a several physical processes which are depen- 
dent on conditions of pumping. Under a low excitation 
level, when the one-photon transitions take place, the en- 
ergy relaxation and recombination of photoexcited carri- 
ers were studied with the use of the time-resolved pump- 
probe measurements, see experimental data and theoret- 
ical discussions in Refs. 1 and 2, respectively. Under an 
extremely high pumping, the multi-quantum transitions, 
which cause the harmonics generation and the hybridiza- 
tion of electron-hole states, take place. This regime is not 
investigated completely for graphene, see general consid- 
eration for bulk materials or quantum wells in Sects. 10 
and 56 of Ref. 3. The nonlinear response is also possible 
within the single-photon approach because the Rabi oscil- 
lations of coherent response should take place in graphene 
if a pulse duration is less 100 fs. [4j Beside of this, the 
nonlinear regime of energy relaxation and recombination 
due to the Pauli blocking effect takes place if the photo- 
generation rate is comparable to the energy relaxation or 
recombination rates. A saturation of transient absorp- 
tion, which have been investigated recently, [!, @ is the 
most important manifistation of such a regime. It is be- 
cause this phenomenon was exploited for realization of 
an ultrafast laser with a graphene saturable absorber in 
the telecommunication spectral region. To the best of 
our knowledge, a complete theoretical treatment of the 
saturation mechanism is not performed yet and an inves- 
tigation of this phenomena is timely now. 

In this paper, we consider the temporal nonlinear evo- 
lution of carriers under photoexcitation by the ultrafast 
pulses in mid-IR and visible (or near-IR) spectral regions. 
For the mid-IR pump, one can neglect the relaxation and 
recombination processes (the quasielastic relaxation due 
to acoustic phonons remains uneffective up to nanosec- 



*Electronic address: ftvasko@yahoo.com 



onds) and the transient distribution of nonequilibrium 
carriers is determined by the broadening of interband 
transitions due to an elastic scattering and by the pa- 
rameters of excitation. For the excitation energies, hQ, 
above the optical phonon energy, Tiui n (77 = T, K la- 
bels the phonon modes correspondent to the intra- and 
intervalley transitions), when a cascade emission of the 
optical phonons dominates in relaxation, the transient 
distribution transforms into a set of peaks. The effec- 
tive electron-hole recombination takes place if the lowest 
peak is placed around the half-energy of optical phonon, 
fkOrf/2. So that, the character of response modifies essen- 
tially if the frequency fl varies over the w^-range. 

The saturation process is described within the frame- 
work of the temporally local approach, when the decoher- 
entization time (which determines the broadening of the 
interband transitions) is shorter in comparision with the 
duration of pumping. Spectral and temporal dependen- 
cies of the photoexcited concentration and the response 
on the probe radiation (transmission and absorption co- 
efficients) are presented. The thresholds for saturation 
of response are estimated to be about 0.2 MW/cm 2 , 60 
MW/cm 2 , and 0.6 GW/cm 2 for HQ -0.12, 0.8, and 1.5 
eV, respectively (mid-IR, near-IR, and visible spectral 
regions). These results are dependent on the decoher- 
entization and relaxation mechanisms and the electro- 
dynamics conditions. Their are discussed in compari- 
son with the experimental data for the near-IR pumping 
case ||| . Conditions for the transient negative absorp- 
tion of a probe radiation during cascade emission of opti- 
cal phonons are also analyzed (this phenomenon under a 
steady-state pumping was considered recently [7] in con- 
nection with a possibility of the THz lasing effect). 

The consideration below is organized as follows. The 
temporally local approach for description of the response 
of photoexcited carriers is developed in Sec. II. Spectral 
and temporal dependencies of the response are described 
in Sects. Ill and IV for the cases of excitation in mid- 
IR and visible (or near-IR) spectral regions, respectively. 
A discussion of experimental data, the list of assump- 
tions used, and concluding remarks are given in the last 
section. In Appendix we consider the mechanism of sat- 
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uration caused by the collisionless Rabi oscillations. 



II. TEMPORALLY LOCAL APPROACH 

Since the symmetry of the energy spectrum and scat- 
tering processes for electrons and holes in an intrinsic 
graphene, we describe the phenomena under considera- 
tion by the same distribution functions for the both types 
of carriers, f p t- According to Refs. 2b and 4, the kinetic 
equation for f pt takes form: 



%± = Vpt (l-2f pt ) + J{f t \p) 



(1) 



and it should be solved with the initial condition 
fpt^oo — 0. Here J(ft\p) is the collision integral, which 
is described the relaxation and recombination processes, 
and u pt is the interband generation rate due to the in- 
plane electric field Ewt exp(— iClt) + c.c, where E is the 
field strength, fl is the pumping frequency, and Wt is the 
envelope form-factor of pulse with duration 2t p centered 
at t = 0. Supposing that t p exceeds the dephasing time, 
we have used in Eq. (1) the temporally local approach 
with the separated filling factor, (l — 2f pt ), and with the 
rate of photoexcitaion: 



v P t = v R w 2 A 



2vp — hfl 
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Here v = 10 8 cm/s is the velocity of neutrinolike quasi- 
particles, 27 is the broadening of the interband excitation 
described by the phenomenological factor A(z). Below 
we consider the Lorentzian lineshape of photoexcitation, 
when A(z) = [n(l + z 2 )] and the Gaussian temporal 
envelope w T = Y / 2/7rexp[— (t/r p ) 2 ]. 

The solution of Eq. (1) determines both the photoin- 
duced concentration, which given by the standard for- 
mula 



responsible for the nonlinear behavior of the response un- 
der consideration. 

We restrict ourselves by the the geometry of normal 
propagation of radiation. The relative absorption of 
graphene sheet, £ w t, as well as the reflection and trans- 
mission coefficients, R u t and T wt , arc determined through 
(Jut- Since the energy conservation requirement, [§[ 
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we consider below only the absorption and transmission 
coefficients: 
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Here y/e is the refraction index of a thick substrate (for 
SiC>2 substrate y/e ~1.45 and dispersion of e can be ne- 
glected) and we approximately separated the carrier con- 
tributions using the coefficients £ m w Ana/ [y/e(l + y/e)}, 
T m w 4Ve/(l + v/e + na) 2 , and a w 2na/(l + ^ft + ira) 
with a = e 2 /he. Notice, that the negative absorption 
regime £ wt < takes place if f Pat t > 1/2, under the pop- 
ulation inversion condition (see discussion in Sec. IV). 
At oj = f2 these relations describe the propagation of 
(2) pumping pulse with the time-dependent intensity Sw 2 , 
where S is the maximal intensity. Performing the aver- 
aging of (6) over the pulse duration one obtains 
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where we used J dtw 2 /t p = 1 . Below we solve Eq. 

—00 

(1) and analyze the responses (6) and (7) for different 
parameters of pump and probe radiations. 
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and the transient response on a probe radiation of fre- 
quency u) (oc exp(— iu>t), which is described by the 
dynamic conductivity a U f For the collisionless case 
huj/^f 3> 1, when the parametric dependency on time 
takes place, [t| the real part of <7 ut is written as follows: 



Reer^ = — (1 - 2/ PM)t ), 



(4) 



where p^ — huj/v. The imaginary part of a^t is de- 
termined through Recr^t with the use of the dispersion 
relation and one can check that the carrier-induced con- 
tribution to Imcr wt appears to be weak in comparison 
with (4) for the peak-like distributions of carriers con- 
sidered below. Thus, the only filling factor in Recr^i is 



III. MID-IR EXCITATION 

We consider here the mid-IR pumping case when the 
energy relaxation of carriers is ineffective and J(ft\p) in 
Eq. (1) can be neglected. As a result, the solution of the 
problem (1) takes form: 



fpt — 



dt'i 



•pt> exp 



-2 / dru 



pT 



(8) 



Evolution of such a distribution from zero value at 
t -C — t p to the saturated peak with the maximal value 
fmax — 1/2 is shown in Fig. la versus dimensionless time 
and energy at the pumping intensity S =1 MW/cm 2 . 
Temporal dependencies of f pn t at different S are shown in 
Fig. lb. These calculations were performed for hil ~120 
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meV (pumping by C02-laser), the pulse duration 2t p ~1 
ps, and the broadening energy 7 ~6 meV which is in 
agreement with the mobility data for the case of elastic 
scattering. |j[ The temporally-dependent photoinduced 
concentration n t is plotted Fig. lc for the same parame- 
ters. The saturated concentration versus intensity, which 
is attained at t > t p , is presented for 7 =6 and 12 meV 
in Fig. Id. These dependencies can be fitted as 



bS 



1 + s/s n ' 



(9) 



where b ~6 or 12.2 MW _1 [ns is measured in 10 11 cm~ 2 ] 
and S n ~1.76 or 10 MW/cm 2 for and 7 =6 or 12 meV, 
respectively. 
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FIG. 1: (Color online) (a) Photoexcited distribution f pt 
versus energy vp and dimensionless time, t/r p at mid-IR 
pumping level S —1 MW/cm 2 . (b) Temporal evolution of 
fat = f pnt at S =0.1, 0.3, 1, 3, and 6 MW/cm 2 (marked 
as 1-5). (c) Potoinduced concentration versus t/r p for the 
same conditions as in panel (b). (d) Concentration ns at 
t/rp — > 00 versus S for the different broadening energies 7. 
Dotted curves are correspondent to the fit (9). 

The relative absorption and transmission coefficients of 
a probe radiation of frequency uj are determined through 
f Puit according to Eqs. (6). Spectral and temporal de- 
pendencies of are shown in Fig. 2a for the conditions 
used in Fig. la. Since af Pujt "C 1, the peak of relative 
transmission T wt /T m resembles f pn t presented in Fig. la. 
Here T m =~0.95 is the transmission coefficient without 
for non-doped graphene. The temporally-dependent rel- 
ative absorption and transmission at the pumping fre- 
quency fl and at different S are presented in Figs. 2b 
and 2c, respectively. The saturated values of £s/£ m and 
Ts versus intensity are plotted in the upper and lower 
panels of Fig. 2d for the parameters used in Fig. Id 
(£s and T5 have only a weak dependency on 7). These 
curves can be fitted as 



l + S/S 



hS 



where S «0.2 MW/cm 2 and h «0.09 cm 2 /MW. The 
saturation of £s and Ts takes place at lower threshold 
intensities in comparision to ns, c. f. Figs, lc and 
2d, 2e. Thus, for the pumping range >1 MW/cm 2 the 
one-photon absorption is suppressed and a damage of 
graphene by mid-IR radiation with t p <1 ps is not pos- 
sible. 
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(10) 
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FIG. 2: (Color online) (a) Spectral and temporal dependen- 
cies of relative absorption, £ w t at S =1 MW/cm 2 . (b) Tem- 
poral evolution of £nt at S used in Fig. lb (marked), (c) 
The same as in panel (b) for transmission, Tm- (d) Avaraged 
over pulse absorption and transmission coefficients (upper and 
lower panels, respectively) versus S [dashed curves are corre- 
spondent to Eq. (10)]. 



IV. CASCADE EMISSION EFFECT 

In this section we consider the photoexcitation by vis- 
ible and near-IR radiation, when the cascade emission of 
optical phonons should be taken into account in Eq. (1). 
For the temperatures below the optical phonon energies, 
the spontaneous emission processes are only essential and 
the collision integral is given by the finite-difference form 
(see evaluation in Refs. 2b and 10) 

J Cftb) = [ V P+Pn i 1 ~ fp^ fp+Pn* ( n ) 

n 

~ v p-p v (l — fp—ptjt) fpt — l/ p n —pfp n —ptfpt\ ■ 

Here ij = T, K is correspondent to the intra- and in- 
tervalley transitions with the energy transfer, Tiuj^, and 
the momentum transfer, p v = hu v /v. The last contri- 
bution of Eq. (11) is responsible for the recombination 
process while the first and second terms describe the in- 
terband cascade relaxation of carriers. The relaxation 
rates v p and v p are proportional to the density of states, 
v p ~ v p ~ 9(p)v r ip/h 1 where the characteristic velocities 
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vy.k can be estimated crudely as Vk ~ 2 x 10 6 cm/s 
and i; r « 10 6 cm/s. [T^ Thus, the -ftf-mode emission 
gives a dominant contribution to the relaxation process; 
moreover, the only interband recombination is possible 
in the passive region, < vp < Tilok =170 mcV. Below 
we neglect other relaxation processes, so that a peak-like 
transient distribution of carriers takes place due to the 
negligible phonon dispersion and a narrow distribution 
of photoexcited carriers, under the condition 7 <C Took- 
For the sake of simplicity, the cases of effective or sup- 
pressed recombination, when the lower peak in the pas- 
sive region is placed around or outside the energy hcox 
are considered. It is convenient to analyze calculations 
for the near-IR and visible pumping cases separately. 
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FIG. 3: (Color online) (a) Contour plots of photoexcited 
distributions f pt versus energy and time for pulse duration 
2t p =0.6 ps at pumping level 200 MW/cm 2 and HQ. =850 
meV. (b) Transient evolution of concentration n t and and 
populations of peaks around ~43, 213, and 383 meV (marked 
as ni, fi2, and 77,3, respectively), for the same conditions as 
in panel (a), (c) Evolution of nt for pumping levels 5* =100, 
200, and 300 MW/cm 2 (marked as 1, 2, and 3). Solid and 
dashed curves are plotted for HQ. = 850 meV and 765 meV. 



MW/cm 2 . The carrier concentrations over the peaks 1- 
3 and the total consentration n t given by Eq. (3) are 
shown in Fig. 3b for the same parameters as in Fig. la. 
Since the relaxation rate in Eq. (11) is proportional to 
the density of states, v p oc p, the bottleneck effect takes 
place under the transition between the second and third 
peaks and rii exceeds ni t a- The transient evolution of 
concentration for different S is shown in Fig. 3c where 
the maximal concentration exceeds 10 12 cm -2 at t ~ t p 
and S >0.3 GW/cm 2 . During the further evolution, n t 
decays due to the recombination process. The case of 
the suppressed recombination (fofl =765 meV) is differ- 
ent because of, first, the peaks are shifted below (about 
43 meV) and, second, the decreasing of f pt and n t due to 
recombination is absent. The saturated concentrations 
(dashed curves in Fig. 3c) exceed the peak concentra- 
tions (solid curves in Fig. 3c) by factor ~1.3. Note, that 
nt decreases with increasing of 7 at fixed S (not shown 
in Fig. 3). 





S (MW / cm 2 ) S (MW / cm 2 ) 

FIG. 4: (Color online) (a) Transient evolution of relative ab- 
sorption for ftfi = 850 and 765 meV (solid and dotted curves, 
respectively) at S =50, 100 and 200 MW/cm 2 (marked as 1- 
3). (b) The same as in panel (a) for transmission coefficient, 
(c) Avaraged over pulse absorption versus S [dashed curve is 
correspondent to Eq. (9)]. (d) The same as in panel (c) for 
transmission coefficient. 



A. Near-IR pumping 

First, we consider the three-step cascade processes 
under the near-IR pumping with wavelengths around 
~ 1.5 /im and the pulse duration determined by t p = 
0.3 ps. We consider the regimes of the enhanced or sup- 
pressed recombination supposing Tiil — 850 meV or 765 
meV. For this energy region, the broadening of photoex- 
cited peak is taken as 7 ~18 meV, so that ft/7 <§; r p . The 
numerical solution of Eq. (1) with the collision integral 
(11 ) is performed with the use of the temporal iterations 
(ll| at different S. Figure 2a shows the contour plot 
of the three-peak distribution function f pt for the case 
of efficient recombination (hft = 850 meV) at S =200 



Transient evolutions of the absorption and transmis- 
sion coefficients given by Eq. (6) are shown in Figs. 4a 
and 4b at the different pumping frequencies f2 (solid and 
dashed curves) and at different S. For t/r p <0.5, the 
temporal evolution of £nt and Tot do not dependent on 
the character of recombination. For t/r p <1.5 this evo- 
lution is completely different: a quenching of photore- 
sponse or a steady-state contribution take place for the 
effective or suppressed recombination cases. At S >300 
MW/cm 2 and t/r p ^0 one obtains the saturated absorp- 
tion around £nt ~0.1. The negative absorption takes 
place for a probe radiation with Tilo around the first and 
second peaks. It is because f Pujt > 1/2, see Eq. (4) and 
the contour plot in Fig. 3, where the regions of negative 
absorption are separated by the thick (red) curves. Thus, 
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the negative absorption (and a possible stimulated emis- 
sion of mid-IR radiation) is realized at S >f 00 MW/cm 2 
during time intervals t < 5r p . 

The absorption and transition coefficients averaged 
over pulse duration according to Eq. (7) are shown in 
Figs. 4c and 4d. Since the transient response at \t\ < t p 
does not depend on the recombination mechanism (see 
Figs. 3c, 4a, and 4b), the variation of £g and T$ with Ml 
is less than 5%. These dependencies can be fitted by Eq. 
(10) with the characteristic intensity S «60 MW/cm 2 
and the coefficient h «0.3 cm 2 /GW. 
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FIG. 5: (Color online) (a) Transient evolution of concentra- 
tion m and and populations of peaks around ~0.09, 0.26, 0.43, 
0.6, and 0.77 eV (marked as ni-s, respectively) for pulse dura- 
tion 2t p =0.4 ps at pumping level 0.4 GW/cm 2 and hQ, =1.53 
eV. (b) The same as in panel (a) at TlQ =1.615 eV for peak's 
positions ~0.13, 0.3, 0.47, 0.64, and 0.81 eV marked as nis. 
(c) Evolution of nt for pumping levels S =0.2, 0.3, 0.4, 0.6 
and 0.8 GW/cm 2 (marked as 1-5, respectively) for hQ =1.53 
eV. (d) The same as in panel (c) for Ml =1.615 eV. 



B. Visible pumping 

Next, we consider the visible light pumping, with wave- 
lengths around ~ 0.75 p,m, using the pulse duration 
2t p =0.4 ps and the broadening 7 «34 meV (so that 
fi/'f -C T p ). Supposing Ufl = 1.53 and 1.615 eV for the en- 
hanced and suppressed recombination regimes one arrive 
to the distribution function formed during the five-step 
cascade process. Transient evolutions of the concentra- 
tions over the peaks 1-5 and of the total concentration nt 
at S =0.4 GW/cm 2 are shown in Figs. 5a and 5b for the 
cases of enhanced and suppressed recombination, respec- 
tively. Similarly to the near-IR pumping case, the upper 
peak concentrations decrease fast at t > t p and a maxi- 
mal population of the second peak takes place due to the 
bottleneck effect. Once again, at t < t v the shapes of n t 
are the same for the both cases. At t > t p a quenching 



of n t due to recombination takes place in Fig. 5a while 
there is no a decreasing of n t in Fig. 5b. The temporal 
dependencies of concentration for different S are shown 
in Figs. 5c and 5d for the two recombination regimes 
under consideration. The maximal concentration range 
up to 10 13 cm -2 at t ~ r p and S «1 GW/cm 2 . 
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FIG. 6: (Color online) (a) Transient evolution of relative 
absorption for h£l =1.53 and 1.615 eV, (solid and dotted 
curves respectively) at S =0.2, 0.4, 0.6, 0.8 and 1.2 GW/cm 2 
(marked as 1-5). (b) The same as in panel (a) for transmis- 
sion coefficient, (c) Avaraged over pulse absorption versus S 
[dashed curve is correspondent to Eq. (9)]. (d) The same as 
in panel (c) for transmission coefficient. 

The temporal evolution of £nt and To t at frequency 
f2 (solid and dotted curves are correspondent to the two 
recombination cases under consideration) are plotted in 
Figs. 5a and 5b. By analogy with Sect. IVA, the nega- 
tive absorption regime takes place at S >0.3 GW/cm 2 . 
Beside of this, the conditions £ w t < take place around 
the peak positions at uj < f2 (not plotted, see a similar 
behavior in Fig. 3a); for the first and second peaks the 
negative absorption regime is realized up to t ~ 5r p at 
S >0.1 GW/cm 2 . In addition, at t < 0.5r p the response 
does not dependent on recombination and at t > 1.5r p a 
damping or time-independent response is realized for the 
effective or suppressed recombination. 

The averaged according to Eq. (7) absorption and 
transmission coefficients, which do not depend on the re- 
combination mechanism, are plotted in Figs. 6c and 6d. 
Once again, £s and Tg can be fitted by Eqs. (10) with 
the characteristic intensity S «0.56 GW/cm 2 and the co- 
efficient h sa0.03 cm 2 /GW. Since the departure rate from 
the photoexcited peak increases if hfl grows, the charac- 
teristic intensity S is also increased in the visible spectral 
region in comparison with the near-IR pumping case. 

V. DISCUSSION AND CONCLUSIONS 

To summarize, we have developed the nonlinear the- 
ory of transient response of an intrinsic graphene under 
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the ultrafast interband excitation. Within the local time 
approach, the conditions of saturation of absorption were 
found in the mid-IR, near-IR, and visible spectral regions. 
In addition, we have demonstrated a possibility for the 
stimulated mid-IR radiation due to the bottleneck effect 
during the cascade emission of optical phonons. Our con- 
sideration is based on the set of assumptions about re- 
laxation mechanisms. First of all, the phenomenological 
model for the broadening with the characteristic energy 
7 is used for description of the intersubband transitions. 
In Sects. Ill and IV we estimated 7 from the experi- 
mental data for the departure relaxation rates. 0, 0, III 
Secondary, a simplified description of energy relaxation 
is employed. We neglect the Coulomb scattering which 
is not a dominant relaxation channel at t < r p , so that 
the results for £5 and T5 should not be modified essen- 
tially. But a transient distribution at t 3> t p and a con- 
dition for the negative absorption of a probe radiation 
in the mid-IR region can be modified. Also, a possible 
contribution of the substrate vibration [ijj is not taken 
into account. These points require a special considera- 
tion but, anyway, our calculation gives a lower bound of 
S. The other assumptions (parameters for the electron- 
phonon coupling, conditions for the temporally-local ap- 
proach, and description of the interband response) are 
rather standard for the calculations of the optical prop- 
erties and the relaxation phenomena in graphene. In ad- 
dition, an inhomogenity of pumping, which causes the 
lateral diffusion of carriers, [13[ and a heating of phonons 
[Til may be essential; these phenomena requre a special 
treatment, both experimental and theoretical. 

We turn now to discussion of the experimental data 
available for the near-IR spectral region. || Numer- 
ical estimates for the saturation thresholds and for the 
concentrations of the photoexcited carriers are in a qual- 
itative agreement with the consideration performed. But 
an accurate comparison with the results presented is not 
possible for the two reasons. First, the graphene struc- 
ture was embedded into the laser cavity in [g, Q so that 
the electrodynamical conditions (for a propagated, re- 
flected, and absorbed radiation) were different from the 
simple geometry considered here. Second, the multi-layer 
graphene or the graphene flakes were used, while a single- 
layer graphene case was not under a detailed treatment. 
Thus, a special measurements with the use of the sim- 
plest geometry of a well-characterized sample placed over 
a semi-infinite substrate are necessary. 

In closing, we have analyzed theoretically the con- 
ditions for realization of an efficient graphene-based 
saturable absorber and have performed a comparison 
with the experimental data. More extended treatment 
of this phenomena under near-IR pumping, including 
an above-mentioned special measurements, in order to 
improve an efficiency of the graphene based saturable 
absorber in the lasers for telecommunications. An 
additional study in the mid-IR and visible spectral 
regions should be useful for verification of different 
relaxation mechanisms. 



The author would like to thank E. I. Karp for insightful 
comments. 



Appendix: Rabi oscillations regime 

Below we describe the saturation of the averaged ab- 
sorption and transmission coefficients (7) under an ul- 
trafast pumping for the case when the Rabi oscillations 
conditions are satisfied. [1] The collisionless regime of re- 
sponse is described by the S'-dependent contribution to 
the distribution function 



dt' 

1 - 2/pt = cos I \j I —w t 



(A.l) 



Here the characteristic intensity is given by 



St 



An 



(A.2) 



and Sr ~0.6 MW/cm 2 for CO2 pumping with t p ~0.1 
ps. For the near-IR or visible pumping with r p ~30 
fs, one obtains Sr ~0.3 or 1 GW/cm 2 . Notice, that 
Sr oc (fl/rp) 2 and (A.l) is not dependent on any other 
parameter if r p is shorter than the dephasing relaxation 
time. 




Vs/s. 



FIG. 7: (Color online) Normalized absorption coefficient given 
by (A.3) versus dimensionless intensity S/Sr. Dotted curve 
presents a monotonic fit. 

Substituting the distribution (A.l) into Eqs. (6) and 
(7) one obtains the following analytical expression for the 
relative absorption 



sm 



dt 



-w t cos 




(A.3) 



while the transmission coefficient is given by Tg T m (l — 
a£s) with a = \/e(l + y/e)/2. In Fig. 7 we plot the func- 
tion £s/£ m versus dimensionless intensity S/Sr and the 
oscillating character of response at S/Sr >4. The oscil- 
lations appears due to the dynamic inversion of transient 
population, see Ref. 4. The fit of (A.3) at S/S R < 1 



is given by Eq. 
S = 2S R . 



(10) with the characteristic intensity 
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